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Abstract — In this work, the Power Summation Method is 
adapted to be applied in distribution systems that include 
distributed generation. The voltage control at the generation 
buses is treated as an optimization problem, which uses the 
Newton-Raphson method as the minimization technique. The 
algorithm is iterative and includes the reactive and voltage limits 
constraints for each generator. The accuracy and computational 
efficiency of the proposed method were evaluated by comparing 
its performance with the Newton-Raphson conventional power 
flow. The test system was a 36 buses feeder with two generators 
installed. 
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I. INTRODUÇÃO 


ários aspectos de interesse das companhias distribuidoras 

fazem com que a implantação da geração distribuída 
(GD) nos sistemas de distribuição seja cada dia mais comum. 
Considera-se como GD aquela que é levada a efeito com 
unidades geradoras de pequena escala (até 10 MW) dispersas 
ao longo de um sistema de distribuição [1]. Como regra geral, 
aproximar a geração das cargas traz benefícios à distribuição e 
transmissão, pois além de permitir que as empresas 
posterguem investimentos em infra-estrutura, é possível 
reduzir perdas e melhorar o perfil de tensão nos alimentadores. 
Com uma configuração adequada, a GD pode também 
aumentar a confiabilidade do sistema e melhorar a qualidade 
da energia entregue ao consumidor final. 

Apesar dos benefícios inerentes à GD, integrar pontos de 
geração a sistemas que, costumeiramente, possuíam apenas 
uma direção de fluxo de potência não é tarefa fácil. Novos 
problemas, como os conseqüentes distúrbios no sistema de 
proteção e a possível formação de ilhas (regiões isoladas do 
alimentador supridas por uma unidade geradora), são novos 
desafios a serem enfrentados. Outra questão que desperta 
interesse dos pesquisadores é o dimensionamento e o 
posicionamento ótimo das unidades geradoras no alimentador 
[2], assim como a combinação ótima entre energia proveniente 
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da subestação e energia proveniente dos geradores, sob vários 
aspectos, como o econômico [3], de redução de perdas [4] e de 
melhoria da confiabilidade [5]. 

Inicialmente, os métodos de fluxo de carga baseados na 
matriz de impedância nodal, como o Newton-Raphson 
(FCNR) e o desacoplado rápido, foram aplicados aos sistemas 
de distribuição. Entretanto, devido à sua elevada razão R/X e à 
sua estrutura predominantemente radial ou pouco malhada, as 
redes de distribuição se enquadravam na categoria de sistemas 
mal-condicionados para estes métodos matriciais, o que trazia 
dificuldade de convergência em alguns casos. Visando 
contornar estas deficiências, novos métodos de fluxo de carga 
para sistemas de distribuição foram propostos. Entre eles, o 
Método da Soma de Potência (MSP) [6-7] é um dos mais 
difundidos. Este método leva em consideração as 
características do sistema de distribuição para desenvolver um 
algoritmo de fluxo de carga robusto, eficiente e não-matricial, 
alcançando a convergência até em casos em que os métodos 
matriciais divergiriam. Além do MSP, foram desenvolvidos o 
Método da Soma de Corrente (MST) e o Método da Soma de 
Impedância, todos possuindo em comum o fato do algoritmo 
de solução seguir duas direções de cálculo ao longo das barras 
do alimentador, por isso sendo chamados de métodos 
backward/forward. 

Com a difusão da GD, torna-se necessário uma adaptação 
dos métodos backward/forward para a inclusão de barras de 
tensão controlada (PV) nos alimentadores. Em [8], o MSI é 
modificado para possibilitar tanto o fechamento de alguns 
anéis no alimentador, quanto para a inclusão de barras PV. 
Nesse trabalho, como as barras em que há fechamento de anel 
são tratadas de forma diferente das barras PV, o acréscimo 
destas piora a eficiência computacional do algoritmo. O uso do 
MSP em vez do MSI em [9] permitiu simplificar o tratamento 
dado às barras PV, reduzindo consideravelmente o esforço 
computacional do algoritmo. Portanto, houve um avanço do 
método de solução apresentado em [8]. Em [10], foi proposta 
uma nova forma de atualização das correntes nas barras PV, 
visando melhorar a eficiência e precisão dos métodos 
backward/forward. 

O principal objetivo deste trabalho é modificar o Método da 
Soma de Potência para estender sua aplicabilidade a sistemas 
com geração distribuída. O controle da tensão nas barras com 
geração é abordado como um problema de otimização. Nas 


Seções II e III o problema é formulado matematicamente e é 
apresentado o método de solução utilizado. Na Seção IV 
encontram-se os resultados e análises da comparação entre o 
método proposto e o FCNR, aplicados a um alimentador radial 
de 36 barras cujos dados foram obtidos em [11]. Por fim, na 
Seção V são apresentadas as conclusões do trabalho. 


Il. FORMULAÇÃO DO PROBLEMA 


O problema de incluir geração distribuída no método da 
soma de potência tem a seguinte abordagem: calcular as 
potências reativas que devem ser injetadas por todos os 
geradores instalados no sistema, objetivando igualar as tensões 
nas barras de geração aos valores especificados para cada uma 
delas. De forma geral, o que se pretende é determinar o 
mínimo de uma função de erro f que tem como argumento uma 
configuração factível dos geradores. 

Considere um alimentador de n barras, em que m delas são 
de geração e as demais são de carga. Em cada uma das barras 
de geração existe um gerador instalado, pertencente ao 
conjunto ®={G;, Gz..., Gm}. Para todas estas barras de geração, 
há valores de tensão especificados 
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que devem ser mantidos às custas da injeção de potência 
reativa 


Qs =106,.06, 26,1 2) 


m 


por parte dos geradores. Os limites máximo e mínimo de cada 
um deles 
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devem ser respeitados, logo 
Q min S Qg S Q max È (4) 


A tensão em uma barra de geração depende dos reativos 
injetados pelo gerador instalado nesta barra, como também dos 
reativos injetados pelos demais geradores existentes no 
alimentador 

- T 
Va =[VG Ve, eso] e (5) 
Esta tensão deve se localizar entre os valores mínimo e 
máximo de tensão que asseguram a operação segura de cada 
gerador 


V min = [V ming Vmino, o Vming,, I e a 
Vmax = [Vmax º Vmaxo, = Vmaxg,, 1" 
ou seja, 
Vmin < VG $ Vmax - (7) 
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O problema pode, então, ser formulado matematicamente 
da seguinte maneira: 


min f(QG) 
a a E <Q < Q mnax (8) 
sujeito a 
Vmin S Vo = V inak 


NI. MÉTODO DE SOLUÇÃO 


Da maneira como foi formulado o problema, o MSP 
convencional, isto é, sem a inclusão das barras de tensão 
controlada, transforma-se em um bloco funcional do algoritmo 
principal de solução. Para a representação dos trechos do 
alimentador, foi utilizado o mesmo equivalente elétrico de [7] 
para duas barras adjacentes, com apenas uma adaptação: a 
barra terminal de cada trecho passa a possuir não só demanda 
de potência, mas também a opção de geração ativa e reativa, 
conforme apresentado na Fig. 1. 

A solução é encontrada formulando-se um problema de 
otimização iterativo, no qual, a cada iteração, buscam-se 
valores de reativos que são injetados pelos geradores do 
sistema com o objetivo de minimizar uma determinada função. 
Esta, por sua vez, deve ser proporcional aos erros entre a 
tensão medida e a tensão especificada em todas as barras de 
geração. Desta forma, a função objetivo f pode ser definida 
como: 


FQ) =} Ve, Qe) -Vg 1 (9) 
i=l 


Testes prévios de busca exaustiva, realizados em diferentes 
alimentadores de teste e com várias combinações para os 
limites de reativos dos geradores, apontaram que a função f 
apresenta um comportamento quadrático, semelhante ao 
apresentado na Fig. 2 para o caso de dois geradores. O gráfico 
mostra a variação de f quando os reativos injetados pelos 
geradores variam ao longo da faixa de valores possíveis, desde 
Qnin até Qua Tendo em vista este comportamento de f, o 
método de Newton-Raphson (NR) foi o escolhido para a 
procura de seu valor mínimo. 

Sendo b o número de barras PV existentes no sistema, a 
condição necessária, embora não suficiente, para que f seja 
mínima é que 

of 


=0 ,i=12,..,b. 
306, 


(10) 


O método de NR é, então, aplicado à resolução do sistema 
de equações não-lineares proveniente da condição estabelecida 
em (10). 
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Fig. 1. Equivalente elétrico dos trechos do alimentador 
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Fig. 2. Comportamento de f para o caso de dois geradores 
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As derivadas parciais primeiras de f são funções não- 
lineares do vetor das variáveis de controle QG e podem ser 
linearizadas, de acordo com a série de Taylor, da seguinte 
forma: 
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para i=1, 2, ..., b, em que 


AQG = [406,. 406, +- A06, 1 (12) 


é o vetor das atualizações de potência reativa em cada um dos 
geradores das barras PV, variável a ser determinada. 

O sistema de equações formado por (11), em uma notação 
mais compacta, resulta na seguinte equação matricial: 
































d d 9? 
= -- f T AQG. 13) 
Oolosmos alas élo, 
em que 
T 
QG |906, 906, d06, 
e 
92 f 92 f 92 f 
05, 906,906, d06,906, 
2, |--A d° f d? f 
302 | 906; 306, 906,906, (15) 
QG E 2 
M M O M 
3f 3f 3f 
06,906, 906,906, 30% 
lx b 





3 


O problema está então equacionado para ser resolvido 
mediante um processo iterativo de determinação do vetor AQg 
que anule o lado esquerdo de (13). Isso implica, em cada 
iteração, resolver o sistema de equações lineares: 


QG 


df 
0% 





G (16) 
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Como a função f não tem uma forma analítica fechada, 
torna-se necessário aproximar as derivadas parciais 
numericamente. Por definição, os elementos do vetor gradiente 
em (14) são dados por: 


F pp 1 Q0 +10) 
Os, 





(7) 


e 50 e 


sendo ô; um vetor em que o elemento da posição j é igual a £ e 
todos os demais são nulos. Desta forma, tomando-se € finito, 
mas com um valor suficientemente pequeno, a equação (17) 
pode ser aproximada por: 


of 1(06+5;)-f(O6) 
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Para se determinar a matriz Hessiana em (15), calcula-se a 
derivada parcial primeira de (18): 
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Substituindo-se (18) em (19), obtém-se a expressão geral: 
3f 1 
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— f(Qg +8) + f(Qo)l , j,k =1,2,...,b 


ou, para o caso particular em que j=k: 
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De fato, cada cálculo de f em (18), (20) e (21) requer a 
execução do MSP convencional seguidas vezes, com cada 
gerador injetando no sistema os reativos indicados como 
parâmetros da função. A cada execução, os valores das tensões 
nas barras de geração são coletados e utilizados no cálculo de 
(9). A cada iteração, na construção da matriz Hessiana e do 
vetor gradiente em (16), o número máximo de vezes que o 
MSP convencional será executado Nmax é proporcional ao 
número b de barras PV, sendo dado por: 


N max = b(4b+ 1). (22) 


Existem várias formas de se aproximar numericamente a 
derivada de uma função. A aproximação por diferença 
progressiva de (18), embora de precisão menor que a 
aproximação por diferença central, foi escolhida porque exige 
um número menor de execuções do MSP convencional a cada 
iteração. Isto acontece porque, utilizando-se (18), alguns 
cálculos de f para a determinação de um dado elemento do 
sistema podem ser reutilizados na formação de outros 
elementos a cada iteração. Na Fig. 3, pode-se verificar a 
diminuição do número de execuções do MSP convencional ao 
se reaproveitar valores de f. 

O algoritmo de solução completo do problema é 
apresentado no fluxograma da Fig.4. Após receber as 
informações necessárias dos trechos do alimentador e dos 
geradores, todos os geradores passíveis de permanecerem 
ligados no sistema recebem suas estimativas iniciais de 
reativos injetados. Neste trabalho, a estimativa inicial utilizada 
foi a média entre os valores mínimo e máximo de cada 
gerador. 

O processo iterativo inicia com a construção das matrizes de 
derivadas parciais em (16), somente fazendo parte da 
formação destas matrizes, as g barras PV do alimentador. 
Construída a equação matricial de (16), são então obtidas e 
implementadas as atualizações de potência reativa para cada 
uma das barras PV do sistema, considerando-se as restrições 
do problema. 

O número g de barras PV é variável ao longo do processo 
iterativo. Inicialmente, todas as barras de geração (onde 
existem geradores instalados) são barras PV, ou seja, g=m. 
Entretanto, conforme indicado no fluxograma da Fig. 4, a cada 
iteração algumas barras PV podem ser transformadas em 
barras PQ, tal qual acontece no método de fluxo de carga de 
Gauss-Sidel. Isto acontece quando a atualização dos reativos 
em uma dada barra PV ultrapassa os limites mínimo ou 
máximo de reativos do gerador. Esta barra é, então, 
transformada em barra de carga com sua potência reativa 
tornando-se constante e fixada no valor limite que foi 
ultrapassado. 

A transformação entre tipos de barras é reversível, ou seja, 
uma barra PV que foi transformada em barra de carga em uma 
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Fig. 3. Comparativo do número de execuções do MSP convencional 
com e sem reaproveitamento de valores na construção de (16) 





Fig. 4. Fluxograma do algoritmo de solução 


dada iteração pode voltar ao seu estado original em uma 
iteração subsegiiente. Isto acontece sempre que houver a 
possibilidade de variar os reativos injetados pelo gerador desta 
barra no intuito de aproximar a tensão na mesma do seu valor 
especificado. 

Após realizar as transformações necessárias entre tipos de 
barras, é verificada a convergência do algoritmo. Caso a 
diferença entre a potência reativa injetada por todos os 
geradores em duas iterações sucessivas seja menor que a 
tolerância especificada, considera-se que o processo convergiu 
e segue-se para a impressão dos resultados. Caso contrário, 
uma nova iteração do processo iterativo é iniciada com a 
atualização dos parâmetros de (16), considerando para tal a 
nova configuração de barras PV e PQ provenientes dos testes 
realizados na iteração anterior. 


IV. ANÁLISE DOS RESULTADOS 


Com o propósito de validar os resultados e avaliar a 
eficiência computacional do algoritmo proposto na Seção MI, 
o FCNR convencional foi tomado como referência. Os dois 
métodos foram aplicados na solução do fluxo de carga de um 
alimentador teste de 36 barras no qual foram instalados dois 
geradores. Este sistema é originário do alimentador teste de 37 
barras do IEEE em [11], agrupando duas das barras existentes 
e multiplicando-se todas as cargas por um fator igual a quatro. 
Este acréscimo nas cargas teve como objetivo aumentar o 
carregamento do alimentador e justificar a inclusão dos dois 
geradores. O sistema possui apenas uma subestação de 4,8 kV 
e os geradores foram instalados nas barras 7 e 19. Na Fig. 5 é 
mostrado o diagrama esquemático do alimentador; já na 
Tabela I encontram-se listados os dados de cada um dos 
trechos; por fim, na Tabela II são apresentados os dados dos 
dois geradores. 


Fig. 5. Diagrama esquemático do alimentador teste de 36 barras 


Aplicando-se o algoritmo proposto ao alimentador teste 
foram obtidos os mesmos resultados que os apresentados pelo 
FCNR. Os resultados do fluxo de carga para cada um dos 
trechos encontram-se também na Tabela I, enquanto que os 
referentes aos geradores estão apresentados na Tabela III. 

Para avaliar a eficiência computacional dos métodos, os 
dois algoritmos foram executados várias vezes em um 
microcomputador de 2 GHz de freqgiiência de relógio e 480 
MB de memória RAM. O tempo médio de execução do FCNR 
foi de 0,6251 s, enquanto que para o algoritmo proposto este 
tempo foi de 0,4079 s, ou seja, em média 34,75% menor. Em 
alimentadores com um maior número de barras essa diferença 
tende a aumentar, pois o FCNR necessitará operar com 
matrizes cada vez maiores. 

Além do menor tempo de processamento apresentado pelo 
algoritmo proposto, este também apresenta a vantagem de ser 
mais robusto para o cálculo do fluxo de carga em 
alimentadores que apresentam grande carregamento. Essa 
vantagem é inerente aos métodos backward/forward, em 
comparação aos métodos matriciais como o FCNR. 

Analisando os resultados provenientes do fluxo de carga, 
observa-se que as perdas ativas no alimentador foram 
reduzidas com a participação dos dois geradores. Antes elas 
eram de 71,73 kW, passando para 35,50 kW após a instalação. 
Os geradores também ajudaram a melhorar o perfil de tensão 
do alimentador. A menor tensão encontrada do alimentador 
sem barras de geração era de 4,53 kV (0,94 p.u.) e passou a 
ser de 4,61 kV (0,96 p.u.) após a instalação. 


V. CONCLUSÕES 


Neste trabalho, o MSP foi adaptado para aplicação em 
sistemas de distribuição com geração distribuída. O problema 
de manter as tensões especificadas nas barras de geração foi 
formulado como um problema de minimizar uma função erro 
de característica quadrática. Para isto, foi utilizado o método 
de Newton-Raphson, em que os elementos do vetor gradiente e 
da matriz Hessiana são obtidos com aproximação das 
diferenciais por diferenças finitas. 

O algoritmo proposto é iterativo e inclui as restrições de 
injeção de reativos e de limites de tensão de operação dos 
geradores. A convergência em condições de extremo 
carregamento, que é característica no MSP original, foi 
mantida e, embora matrizes e derivadas numéricas sejam 


TABELA I 
DADOS DE ENTRADA E DE SAÍDA DO ALIMENTADOR 












































De Para R(Q) X(O) P(MW) O(Mvar) VkV) (°) P(MW) Q(Mvar) 
O 1 0,0368 0,1404 0,8400 0,4200 4,7589 -0,5662 1,8446 0,9033 
1 2 0,0375 0,0805 0,0000 0,0000 4,7428 -0,7256 1,0025 0,4789 
2 3 0,0515 0,1107 0,0000 0,0000 4,7287 -0,8764 0,6682 0,2906 
3 4 0,1017 0,0728 0,1133 0,0533 4,6949 -0,9461 1,1582 0,5637 
4 5 0,0339 0,0243 0,0000 0,0000 4,6847 -0,9671 1,0428 0,5089 
2 6 0,0610 0,0437 0,1133 0,0533 4,7409 -0,7248 0,0952 0,0734 
6 7 0,0881 0,0631 0,0000 0,0000 4,7410 -0,7173 -0,0181 0,0201 
7 8 0,1356 0,0971 0,1133 0,0533 4,7238 -0,7519 0,4413 0,2165 
2 9 0,1167 0,0435 0,0000 0,0000 4,7359 -0,7186 0,2374 0,1120 
9 10 0,0933 0,0348 0,1240 0,0587 4,7330 -0,7156 0,1240 0,0587 
9 11 0,0700 0,0261 0,1133 0,0533 4,7339 -0,7166 0,1133 0,0533 
7 12 0,0233 0,0087 0,0507 0,0240 4,7400 -0,7164 0,1641 0,0774 
12 13 0,1517 0,0565 0,1133 0,0533 4,7357 -0,7121 0,1133 0,0533 
8 14 0,2683 0,1000 0,0000 0,0000 4,7055 -0,7285 0,2708 0,1347 
14 15 0,2217 0,0826 0,0560 0,0280 4,7024 -0,7244 0,0560 0,0280 
8 16 0,1017 0,0728 0,0000 0,0000 4,7222 -0,7550 0,0560 0,0280 
16 17 0,0817 0,0304 0,0560 0,0280 4,7211 -0,7535 0,0560 0,0280 
14 18 0,0350 0,0130 0,2147 0,1067 4,7036 -0,7261 0,2147 0,1067 
3 19 0,0700 0,0261 0,0560 0,0280 4,7376 -0,8930 -0,4987 -0,2790 
19 20 0,0475 0,0340 0,0560 0,0280 4,7338 -0,9004 0,2801 0,1400 
20 21 0,0817 0,0304 0,0560 0,0280 4,7326 -0,8989 0,0560 0,0280 
20 22 0,0583 0,0217 0,1680 0,0840 4,7313 -0,8972 0,1680 0,0840 
5 23 0,1017 0,0728 0,1133 0,0533 4,6814 -0,9745 0,1133 0,0533 
5 24 0,0542 0,0388 0,0000 0,0000 4,6701 -0,9970 0,9268 0,4537 
24 25 0,0933 0,0348 0,0560 0,0280 4,6688 -0,9952 0,0560 0,0280 
24 26 0,0542 0,0388 0,1133 0,0533 4,6565 -1,0253 0,8685 0,4240 
26 27 0,0949 0,0680 0,0560 0,0280 4,6357 -1,0682 0,7520 0,3684 
27 28 0,1085 0,0777 0,1867 0,0933 4,6190 -1,1024 0,5246 0,2577 
28 29 0,0678 0,0486 0,1680 0,0827 4,6123 -1,1166 0,3375 0,1641 
29 30 0,0678 0,0486 0,0000 0,0000 4,6090 -1,1239 0,1694 0,0813 
30 31 0,0678 0,0486 0,0560 0,0280 4,6079 -1,1261 0,0560 0,0280 
27 32 0,1517 0,0565 0,0000 0,0000 4,6291 -1,0608 0,1694 0,0814 
32 33 0,3733 0,1391 0,0560 0,0280 4,6238 -1,0536 0,0560 0,0280 
32 34 0,0583 0,0217 0,1133 0,0533 4,6275 -1,0590 0,1133 0,0533 
30 35 0,0583 0,0217 0,1133 0,0533 4,6073 -1,1221 0,1133 0,0533 
TABELA II 
DADOS DOS GERADORES 

Barra (G) Pa(MW) Jomin fPmax Vmin(pu) V (pu) Vmax (pu) 
7 0,625 0,85 0,95 0,8500 0,9877 1,1500 

19 0,835 0,82 0,90 0,8500 0,9870 1,1500 

TABELA M 
RESULTADOS DOS GERADORES 

Barra (G;) Qmin (Mvar) Q (Mvar) Qmax (Mvar) V (pu) V (pu) 
7 0,205428 0,274793 0,387340 0,9877 0,9877 

19 0,404409 0,447209 0,582834 0,9870 0,9870 








utilizadas, não houve grande impacto na complexidade do 
algoritmo. O método foi utilizado na resolução do fluxo de 
carga de um sistema de 36 barras e seus resultados 
comparados com os do FCNR. As comparações confirmaram a 


exatidão do algoritmo proposto e também sua maior eficiência 
computacional. 

A principal deficiência observada na estratégia proposta é a 
necessidade de aproximar numericamente várias derivadas. 
Para atenuar essa característica, estão sendo investigadas 
formas de melhorar o procedimento de construção do sistema 
de equações do FCNR, assim como outras estratégias de 
minimização que simultaneamente evitem a construção de 
derivadas numéricas e incorporem as restrições do problema. 
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